Renormalized field theory of collapsing directed randomly branched polymers 
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We present a dynamical field theory for directed randomly branched polymers and in particular 
their collapse transition. We develop a phenomenological model in the form of a stochastic response 
functional that allows us to address several interesting problems such as the scaling behavior of the 
swollen phase and the collapse transition. For the swollen phase, we find that by choosing model 
parameters appropriately, our stochastic functional reduces to the one describing the relaxation 
dynamics near the Yang-Lee singularity edge. This corroborates that the scaling behavior of swollen 
branched polymers is governed by the Yang-Lee universality class as has been known for a long time. 
The main focus of our paper lies on the collapse transition of directed branched polymers. We show 
to arbitrary order in renormalized perturbation theory with e-expansion that this transition belongs 
to the same universality class as directed percolation. 

PACS numbers: 64.60.ae, 05.40.-a, 64.60.Ht, 64.60. Kw 
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I. INTRODUCTION 

It is well known that a single linear (non-branched) 
polymer consisting of N monomers in a diluted solution 
with a good solvent where it is in a swollen conformation 
with gyration (Flory) radius R N ~ N USAW , v SAW > 1/2, 
undergoes a phase transition to a collapsed "globule" 
with Flory radius Rn ~ iV 1 ^ (where d is the dimen- 
sionality of space) if the solvent quality deteriorates as it 
usually does when the temperature of the solution is low- 
ered below the so-called 0-point. Many theoretical tools 
have been successfully applied to linear polymers, and 
considerable progress has been made in understanding 
this collapse transition [l[ . 

In contrast, much less is known about collapse tran- 
sition of randomly branched polymers, which are usu- 
ally assumed being macroscopically either isotropic or 
directed. As far as the isotropic case is concerned, a num- 
ber of numerical studies has been performed over the last 
20 years or so d d 0, H !, 0] . The picture that arises 
from these studies is much more complex than that for 
the collapse of linear polymers. The basic reason for this 
added complexity is the possibility for introducing more 
than one fugacity to drive the collapse. It appears that 
the line of collapse transitions in the phase diagram con- 
sist of 2 qualitatively different parts. One part of the line, 
called the line of ^-transitions, corresponds to continuous 
transitions with universal critical exponents from swollen 
polymer configurations with tree-like character to com- 
pact coil-like configurations. The other part of the tran- 
sition line, the line of ^'-transitions, corresponds to the 
collapse of swollen foam- or sponge-like polymers with 
many cycles of bonds between the monomers to vesicle- 
like compact structures. If this transition is assumed to 
be continuous, one finds nonuniversal exponents (2j. The 
two different parts of the collapse-transition line are sep- 
arated by a higher multicritical point which belongs to 



the isotropic percolation universality class. One of the 
open questions regarding the phase diagram is the ex- 
istence of a possible further transition line between the 
configurations of the collapsed polymers. 

As far as directed branched polymers (DBP) and 
their collapse (CDBP) is concerned, the theoretical pic- 
ture [J, H, 0, [l(| El is somewhat clearer than in the 
isotropic case, in particular in 1 + 1-dimensions (one 
transversal and one longitudinal; in the following, we will 
denote the full space dimension as D = d + 1). Noting 
that branched polymers, lattice animals, and lattice trees 
belong to the same universality class, one can use the re- 
sults of Dhar on collapsing directed strongly embed- 
ded or site animals. This way, Henkel and Seno Q have 
shown that there is only one type of ^-transitions which 
describes the collapse of directed branched polymers in 
(D = 2)-dimensions. This collapse transition belongs 
to the directed percolation (DP) universality class. Fur- 
thermore, Dhar has shown in Ref. that the collapse 
transition is also of directed percolation type in D = 3 
if a special relation between the potentials (fugacities) 
of directed site animals holds. Flory theory was applied 
by Redner and Coniglio @. They correctly derived the 
upper critical dimension d c = 4, but they found criti- 
cal exponents for the CDBP transition that are different 
from the DP exponents which, however, is incorrect as 
we shall see. 

Whereas a comprehensive field theory for swollen (ex- 
tended) polymers has been existing for quite some time, 
the field theory for the collapse of branched polymers 
is much less developed. As far as we know, there ex- 
ist to date no field theory for the directed case. For the 
isotropic case, there is the seminal work by Lubensky and 
Isaacson [l2j and Harris and Lubensky However, 
it turns out that these papers, as far as they consider 
the collapse transition, contain a fundamental error in 
the renormalization procedure (they overlook a required 
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renormalization) , and as a consequence the long-standing 
1-loop results for the collapse transition are not entirely 
correct flij . In an upshot, one may say that there are 
important open questions and unresolved issues in the 
theory for the collapse of directed and isotropic branched 
polymers. 

Here and in a following publication (3], we develop 
dynamical field theories for the collapse of directed and 
isotropic branched polymers, respectively. The idea be- 
hind these theories is to start out with stochastic activa- 
tion processes which lead to percolation clusters and then 
to exploit the well know connection between branched 
polymers and lattice animals, i.e., cluster of a given mass. 
We concentrate on large animals below the percolation 
point, and we equip the theory with enough parameters 
to allow for triciticality which corresponds to the collapse 
of large branched polymers. As alluded to above, the 
directed case has the benefit of being simpler than the 
isotropic case. It allows us to learn about the fundamen- 
tal structure of dynamical field theories for the collapse of 
branched polymers and thereby to sharpen our tools for 
the more complicated isotropic case. The most important 
concrete result of the present paper is that we show to 
arbitrary order in renormalized perturbation theory with 
e-expansion that the CDBP transition is of DP-type. 

The outline of the remainder of this paper is as fol- 
lows: In Sec. HIl we develop our field theoretic model. 
In Sec. IIIBI we analyze our model in mean-field the- 
ory. We determine the mean-field phase diagram of DBP 
and discuss their scaling behavior at the collapse transi- 
tion ignoring fluctuations. In Sec. IIII1 we briefly discuss 
swollen DBP by making contact to established theories. 
In Sec. HV[ we present the core of our field theoretic anal- 
ysis of the collapse transition. We discuss the scaling 
invariances of our model at this transition and their con- 
sequences. We calculate the counter-terms required to 
renormalize our theory in a 1-loop calculation and by 
using Ward identities. We set up and solve renormal- 
ization group equations that provide us with the scaling 
form of the equation of state and correlation lengths. In 
Sec. |V1 we translate our field theoretic results for the 
collapse transition into a form that is more commonly 
used in polymer theory, and we compare them to numer- 
ical simulations. In Sec. IVI1 we give a few concluding 
remarks. There are tow appendixes. In Appendix [Al 
we present some details of our 1-loop calculation, and in 
Appendix [Bl we discuss essential singularities. 



II. TOWARDS A DYNAMICAL FIELD THEORY 
OF DIRECTED BRANCHED POLYMERS 

A. A generalized directed percolation process 

In this subsection, we develop a field theoretic stochas- 
tic functional [U, El ES El for DBP based on very 
general arguments alluding to the universal properties of 
a corresponding Markoffian stochastic percolation pro- 



cess. We consider polymers as clusters in d transversal 
and 1 longitudinal (time) directions. Below the perco- 
lation threshold all of the clusters generated by such a 
process are finite. The distribution function for the large 
clusters decays exponentially with their mass (number 
of monomers). It is well known fl9| that untypical very 
large clusters with linear sizes essentially larger than the 
correlation lengths are generated as rare events. The sta- 
tistical properties of these fractal clusters belong to the 
universality class of lattice animals. Thus, they also form 
the prototype model for single randomly branched poly- 
mers in a dilute solution. Allowing for a further mech- 
anism in the percolation process which favors contacts 
of the monomers, we introduce a possible tricritical in- 
stability that leads to the collapse of the large clusters 
describing the branched polymers. 

A Langevin equation which describes directed perco- 
lation processes allowing for tricritical behavior has been 
well known for many years US, HH, [H, [H, [H, [H, 
[2(| [2?| . It is based on the fundamental phenomenologi- 
cal principles of absorbing processes in conjunction with 
a density and gradient expansion. The basic variable or 
field is the density n(x, t) of agents (infected individuals), 
which models the fractal monomer-density of the poly- 
mer generated from a given source h(x, t). This choice of 
field is possible because the stochastic process generates 
clusters such that the numbers of monomers, branching-, 
and end-points are all of the same order which is in con- 
trast to the case of linear polymers. In other words, the 
density is a proper field for the current problem because 
branched polymers are really "fur-bearing animals" . The 
Langevin equation in Ito-interpretation is given by 

A _1 n = R(n) n + D 1 (n)V 2 n + D 2 {n) (Vn) 2 + ... + h + (, 

(2.1a) 

where A is a kinetic coefficient, and the noise correlation 
reads 



C(x,t)C(x',t') = A- 



Q(n)n(x,t)+... J(x-x') 6(t-t') . 

(2.1b) 

We assume Gaussian noise which we can for our purposes 
without loss of generality because higher order noise does 
not change the universal properties of the process. The 
leading terms of the expansions are Di(n) = l+cin+. . ., 
D 2 (n) = -c 2 + .. ., R{n) = -r- g'n/2- fn 2 /6 + . . ., and 
Q(n) = g + . . .. Further terms in the expansions are pos- 
sible but they turn out being irrelevant. In mean- field 
theory, the percolation transition occurs at r = 0. As 
announced above, we restrict oursclf in the following to 
the region r > in which only non-percolating directed 
branched clusters with a typical linear size £ ~ r" 1 / 2 
are generated from a time- and space-localized source. 
However, our primary interest is in the rare events where 
the fractal clusters of linear size essentially larger than 
r~ x / 2 are generated. As long as g' > 0, the second order 
term f'n of the rate R is irrelevant. We permit both 
signs of g' so that our model accounts for a tricritical 
instability. Consequently, we need the second order term 
with /' > for stabilization purposes, i.e., to limit the 
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density n to finite values. The parameters Cj (which are 
irrelevant for the original percolation problem) are as- 
sumed to be positive to smooth density fluctuations. An 
essential property of percolation processes is the linear 
dependence of the noise correlation (|2.1b[) on the density 
if n — ► for g > 0. This property guarantees that the 
percolation process is really absorbing [25| . and guaran- 
tees the description by a density variable. Last but not 
least, all the terms lead to h = if n = and ft = on 
grounds of the absorbing state. 

The approach that we are taking focuses on general 
principles for processes belonging to the same univer- 
sality class and is therefore necessarily phenomenolog- 
ical [HI]. We devise our field theoretic model repre- 
senting the universality class using a purely mesoscopic 
stochastic formulation based on the correct order param- 
eters identified through physical insight in the nature 
of the critical phenomenon. Hence, the stochastic re- 
sponse functional that we are about to derive stays in 
full analogy to the Landau- Ginzburg- Wilson functional 
and provides a reliable starting point of the field theo- 
retic method. 

Alternatively, one might use the so-called "exact" ap- 
proach which, as a central step, consists of reformulating 
a microscopic master-equation for chemical reactions as 
a bosonic field theory on a lattice. For a recent excellent 
review article on this method see (28|. We choose not to 
use this method as our main approach because we feel 
that it treats universal properties not as transparently 
as the purely mesoscopic formulation does. Nonetheless, 
we think that the lattice reactions which would be the 
starting point for the "exact" approach to the current 
problem have some pedagogical value as they can help to 
nurture ones intuition about the epidemic process. Here, 
these reactions comprise the well known reactions leading 
to DP: 

A(x) A(x)+A(x + <5), (2.2a) 

2X{x) X{x), (2.2b) 

X(x) 0, (2.2c) 

where x denotes a lattice point, S denotes a vector to a 
neighboring point, and X(x) stands for an agent (parti- 
cle) at x. Moreover, because we are interested in tricrit- 
icality, we need additional reactions that produce more 
compact clusters: 

X(x-S) + X{x + S) X(x-S)+X(x) 

+ X{x + S), (2.2d) 

3X(x) X(x). (2.2c) 

For illustrations of the reactions defined in Eq. (|2.2p , see 
Fig. [T] and Refs [H IH- Note t hat the single site back- 
reactions, Eqs. (|2.2b[) and (|2.2e[) act against multiple oc- 
cupations of lattice sites and are thus reminiscent of an 
cxcludcd-volume interaction. As mentioned above, we 
will not use the "exact" approach in this paper. However, 
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FIG. 1: Reactions leading to tricritical DP. Open circles indi- 
cate lattice sites. Solid dots symbolize agents occupying these 
sites. 



for the adherents of the "exact" approach, we mention 
without presenting further details that we have verified 
that he bosonic representation of the reaction equations 
(|2.2p leads to the same stochastic functional (called an 
action in that approach) as the purely mesoscopic ap- 
proach. 

Returning to the latter approach, we now recast the 
Langevin equations (|2.ip as a stochastic response func- 
tional [IE EE E3, El 

Jdp = J d d xdt \n{\- l dtn + (r - V 2 )n + 3 -n 2 

+ t-n 3 - Cl nV 2 n + c 2 (Vn) 2 - |fm} . (2.3) 

Here, we have neglected further higher order terms that 
will not result in any relevant contributions to our final 
response functional. Jdp describes the statistics of DP 
clusters generated by the stochastic process (|2.1[) . Hav- 
ing this functional, one can calculate the average of any 
observable that is polynomial in n and n, say h], by 
functional integration with weight exp(— JTdp), 

(O)dp = J P[n,n]C[n,n]cxp(-JDp) 

=:Tr(0[n,n]exp(-J DP )) . (2.4) 
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Here and in the following, path integrations are always 
interpreted in the sense of the (mathematical save) Ito- 
discretization. One has to be careful to set the range 
of the variables of the functional integration correctly as 
there are some subtleties involved. Some authors [28| 
use functional integrals of the type featured in Eq. (|2.4j) 
purely in a restricted sense as a device for generating a 
perturbation expansion. This precautionary view arises 
from a narrow interpretation of the bosonic functional 
integral after elimination of irrelevant terms including 
higher order monomials of the fields which are in gen- 
eral required to ensure the convergence when fields be- 
come large. However these convergence problems can and 
should be avoided by choosing the support for the func- 
tional integration properly, as was shown e.g. by Ciafaloni 
and Onofri [30] in the case of Reggeon field theory (RFT 
which is equivalent to DP). This type of consideration 
leads here to the rule that n is to be integrated along 
the positive real axis whereas the integration of n is per- 
formed along the full imaginary axis. Of course, devia- 
tions in finite regions, as e.g. suggested by saddle-points, 
are possible. 

One of the most important observables in the current 
problem is the mass or total number of monomers of a 
cluster, 



M = X d d xdt n(x, t) , 



(2.5) 



where we have included the kinetic coefficient A into the 
definition for later convenience (if we were not including 
A in the definition of Ai, we had to multiply A4 in sub- 
sequent formulas by A on dimensional grounds). Because 
we are mainly interested in the scaling behavior of a sin- 
gle large polymer with N 3> 1 monomers, we will focus 
in the following on averages that are restricted to clus- 
ters of a given mass N, and we assume without loss of 
generality that this cluster emanates from a weak source 
h(x.,t) = qS(x)S(t) of agents at the origin x = at time 
t = 0. For our general observable 0[n], this leads to 

(0) N T(N) = (S(N - M)0[n] exp(gn(0, 0))) D p 

~ q(5(N - M)O[n]n(0, 0))) D p , (2.6) 

where V(N) is the probability distribution for finding 
a cluster of given mass N with the last equation hold- 
ing only asymptotically for large N and small q. Note 
that the zeroth-order term in the Taylor expansion lead- 
ing to the second line of Eq. (|2.6p vanishes because 
(JF[?i])dp = J-[0] for any functional T of n because DP 
is an absorbing process. Note also, that the formalism 
allows with ease to ask for the probability of clusters 
generated by several sources at different points (ri, ti) by 
inserting more fields n(rj,i,-) in the exponential of the 
average (|2.6[) . Equation (|2.6p implies that the probabil- 
ity distribution for finding a cluster of mass N is given 
by HI 



V(N) = (S(N - M) cxp(qh(0, 0))) Dp 
^q(S(N-M)n(O 1 0)) Bp . 



Note that V{N) is correctly normalized, J dN V(N) = 
(exp((7n))DP = 1, because (.F[7i])dp = .F[0] for any func- 
tional T of n due to causality. 

V(N) is asymptotically proportional (up to nonuni- 
versal amplitudes and an exponential factor /j,q , where 
fiQ is an effective coordination number of the lattice) to 
the lattice animal number An which plays an important 
role in percolation theory. An measures the number of 
directed clusters of size N, weighted by fugacities for dif- 
ferent cluster properties as contacts, loops, and so on. In 
terms of the number of configurations C(N, b, c) consist- 
ing of N sites, b bonds and c not bounded contacts, An 
can be expressed as [l3| 



An = ^C(iV,&,c)A^ 

b,c 



(2.8) 



The fugacities Ai and A2 correspond to the parameters 
of the stochastic functional J7dp- The collapse transition 
occurs for large N when A2 becomes critical because the 
latter rewards contact-rich animals. If we assume univer- 
sality, the probability distribution V{N) and the cluster 
number should be related for N 3> 1 via 



V{N) 



An Ho 



-N 



J2n> An'Hq 



N' 



and conversely 



An vtfV(N) 



(2.9) 



(2.10) 



a relation that we shall use later on. 

Now we turn to the signature of the collapse transition. 
Chemically, in a diluted solution the transition occurs 
when the hard core repulsion of the monomers is exactly 
compensated by their weak effective attraction, i.e., when 
the second Virial coefficient A^ 2 ' in the expansion of the 
osmotic pressure in powers of the polymer density van- 
ishes 0. Up to a minus sign this second virial-coefficicnt 
is proportional to the space integral of the correlation 
function between two polymers in the solution. Hence, 
in our formalism, the signature of the collapse-transition 
is the vanishing of the space-time integral of the two- 
polymer correlation 



C (2) (TV; x, t) = V {2) (N; x, t) - T {2) (N; 00) , (2.11) 



where 



V {2) (N; x, t) ~(8(N- M) n(x, i)n(0, 0)) 



DP 



(2.12) 



is asymptotically the probability distribution of two large 
clusters of total large mass N generated by two roots at 
(0,0) and (x,i). If the distance of the point (x, t) from 
the origin is substantial bigger than their linear sizes, the 
two clusters decouple, and we have 



(2.7) 



V {2) {N;x,t) ~ V {2) {N;oo) = / dN' V{N')V{N—N') , 

(2.13) 
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i.e., C@'(N; x, t) becomes a cumulant. The second virial- 
cocmcient is proportional to the total measure of this 
correlation times the factor TV 1 / 2 , which comes into play 
if one considers only clusters of the same mass N, 



~ ~N 1/2 (V{N)) A / d d xdtC {2) (N;x,t)q 

~ -N^ 2 X f d d xdt(n{x,t)n(O,0)) { ™ m) . (2.14) 



Instead of working with quantities discussed above as 
functions of TV, it is often more convenient to work with 
their Laplace-transformed counterparts. For V(N) wc 
then have the representation by an inverse Laplace trans- 
formation 



V{N) 



a— too 
cr+ioo 



dz 
2tTi 



c zN (exp(-zM+gn(0,0))) D p 



/ — exp(zA + q$(z) + 0(q 2 )) , (2.15) 

Ja-ioo 27TI 



with $>(z) = (n) z , where (• • • ) z denotes averages taken 
with respect to the new stochastic functional 



Jz = Jdp + zM 



(2.16) 



Note that the normalization condition for V(N) results 
in $(0) = 0. The function Q(z) corresponds to the gen- 
crating function of the animal numbers An, and plays 
in the following the role of an order parameter. Before 
moving on to mean-field theory, we would like to mention 
that a shifted version of the Laplace variable z will play 
later on a role akin to that of an external field. 



B. Mean- field theory 

In this subsection, we discuss various aspects of DBP 
in mean-field theory Our motivation for presenting this 
mean-field theory is twofold, (i) The field theory that 
we are going to present in the following sections is fairly 
involved. Mean-field theory, on the other hand, allows 
us to obtain results with relative ease and may serve the 
reader as a warm up for the more difficult sections to 
come, (ii) At mean-field level, various results for DBP 
are well known. By reproducing these results, we will find 
that our model satisfies important consistency checks. 

Interested in a mean-field approximation in the spirit 
of Landau theory as the starting point of the systematic 
perturbation expansion, wc look for stationary and ho- 
mogeneous saddle-points of J z , which are determined by 



x -i^J± =rh + g'nn+ —hn 2 - -n 2 + z = 0, (2.17a) 
on 2 2 



A _i ^^ =rn + —n 2 + —n 3 — gnn = . (2.17b) 



. t 5Jz 
8h 



These equations are solved by 



\SP 



0, 



$(z) 



n\ sp = (r-r(z))/g 



-(z) = V / 2gz + r 2 . 



(2.18a) 
(2.18b) 

(2.18c) 



The correct sign of the square root is determined from 
the normalization condition $(0) = 0. The equation of 
state for shows a branch-point on the negative real 
axis at z c = —r 2 /2g where r = becomes critical. Note 
that the presence of n\ SP in Eq. (|2.17b|) leads effectively 
to the replacement of the under-critical parameter r by 
the critical parameter r. Hence, this equation leads to 
the usual mean-field equation of tricritical DP [3, [22|, 
[H, [H, |27| . It follows that the saddle-point solution, 
Eqs. (|2.18|) is stable as long as g 1 > 0. If g' becomes 
negative, Eq. (|2.17b|) developes continuously a non-zero 
positive solution n\ SP = —3g'/f at the critical value 
t = 0, signalling the collapse to a compact cluster as a 
continuous phase transition. 

Qualitatively, it is the branch-point singularity that de- 
termines the form of V(N), and higher orders in pertur- 
bation theory, though leading to quantitative improve- 
ment, do not result in qualitative changes of that from. 
To get the asymptotic expansion of "P(-/V), we deform 
the contour of the complex integral in Eq. (|2.15[) as il- 
lustrated in Fig. [2 If Disc $ denotes the discontinuity of 
the function (f> at the branch cut, we get 

P(N) * r dx mScHZc ~ x) e- xN 
Jo 2m 
qe*cN+ q t>(zc) roc Disc$(0 c -x//V) _ T 
dx ; — -e . 



N 



2tti 



(2.19) 



The non-universal factor ge 2cAr+9 *^ 2c - ) is common to all 
the quantities described by Eq. (|2.6p and cancels there- 
fore in all averages (O) n- Note that for large N only the 
small neighborhood of z c matters. In mean-field theory, 
it is easy to calculate V(N) without further approxima- 
tion. We obtain 



V{N) 



-■ill 



exp 



(rN - q)' 



2gN 



qN- 3/2 cxp[-Nr 2 /2g] 



(2.20) 



This distribution has its maximum near the characteris- 
tic value N = N = q/r. As stated above, we are, how- 
ever, mainly interested in the rare events with N 3> Nq 
where the asymptotic forms are valid. Note that the 
second line which gives the large- TV behavior is in per- 
fect agreement with the well known asymptotic result 
7?M)(/V) ~ N- 0{m!> exp(-/Vy/V^ nf) ) with the so-called 
cntropic exponent given in mean-field approximation by 
g(mf) = 3/2 and a non-universal number N^ 1 ^ = 2g/r 2 . 
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Im z 



Re z 



FIG. 2: Deformation of the contour of integration for the 
inverse Laplace transformation. Originally, the integration is 
along the imaginary axis. After deformation, it nestles to and 
leads around the branch cut that terminates at the branch 
point z c on the negative real axis. 



Gaussian fluctuations are governed by the second vari- 
ations of the stochastic functional J z . After Fourier 
transformation in space and time, we have 



ri,i(q, w)mf 
r ,2(q,w)mf 



SJJh, n] 



5nSn 
6J z [h, n] 



Sn5r 



SP 



SP 



^ + A(r + q 2 ), (2.21a) 



A( ff ' + cq 2 )$, (2.21b) 



where c = c\ + cq,. r2,o;mf = — 9 n \sp IS zcr0 - The 
propagator follows from the inverse of Fi x-mf m (q^) - 
representation as 

G?i,i (q, *; z) = 6(t) exp [—A (r(z) + q 2 ) t] , (2.22) 

and the ri-correlation function integrated over space and 
time is 



G ,2(q = 0,w = 0;z) 



r , 2 (o,o) 



mf 



g'(r - t(z)) 



|ri,i(0,0) mf | 2 A. 9 r(z)2 
(2.23) 

The second virial coefficient follows from Eq. (|2.14[) to 
leading order as 



~ N 2 g'r/g. 



(2.24) 



It vanishes at the mean-field collapse transition, that is 
for g'r/g = 0. 

Next, we calculate the average monomer density 
(n(x, i))jv emanating from a source at (0,0). For this 
calculation, we recall Eq. (|2.6[) . When specialized to 
0[n,h] = n(x, t), the second line of Eq. (|2.6p is equal 
to the Laplace transform of (n(x, t)n(0, 0)) z the latter 
being the (space-wise) Fourier transform of Gi j i(q, t; z), 
Eq. (|2.22|) . Collecting, we find that (n(q, i))jv can be 
written as 

— G ljl (q,t;z)exp(zA^) , 

(2.25) 



where the integration path is taken appropriately as out- 
lined above for the probability distribution. The mo- 
mentum integration back to the (x, t) representation is 
straightforward, and we obtain 



(n(x,i))jy 



9*t0(t) 



exp 



x 



g(At) s 

2N 



(2.26) 



Form this result, we can read off directly that the lon- 
gitudinal and transversal radii of gyration scale as R» ~ 



N"£*\ and R i 



(mf) 



N^ A ' ± with the well known mean-field 



animal exponents 

(mf) 
"X II 



Jmf) 



(2.27) 



The stochastic functional for directed branched 
polymers 



The response functional J z = J7dp + zM, Eq. (|2.16[) . 
is form-invariant under three continuous transformations. 
Thus, three parameters of the functional are redundant. 
Here, we exploit the symmetries to eliminate two redun- 
dant parameters and to define scale-invariant effective 
couplings. 

The symmetry transformations are: (i) a rescaling of 
the fields 

f)^a _1 n, n-tan, (2.28a) 
(ii) a mixing of the fields 

h^n + (3n, n^n, (2.28b) 

and (iii) a shift of the response field 

n — > h + 7 , ri — > n . (2.28c) 

We eliminate the redundant parameters r and c by spe- 
cializing the latter two, which do not transform the den- 
sity field n, by 



(r - r) 



1 



(2.29) 



This special shift reminds of the saddle-point, Eq. 
(I2.18bp . and eliminates the uncritical parameter r in fa- 
vor of t. After this shift the fields and their expectation 
values are small quantities for small r, and a perturba- 
tion expansion is appropriate. The parameter r is a free 
parameter at this place, and, of course, we can set this 
parameter to zero. However, it is more economical to 
use it in the renormalization program as a small control 
parameter because it attains at least an additive renor- 
malization. The shift alone introduces besides the bilin- 
ear gradient term ~ fiV 2 n a further quadratic gradient 
term ~ nV 2 n that is thereafter eliminated by the spe- 
cial mixing in Eq. (|2.29[) . This term also attains renor- 
malizations, and to cure UV-divcrgcncies a counter term 
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~ nV 2 n must be reintroduced via the renormalization 
scheme, which we will do as we proceed. 

The transformation (|2 . 29[) suggests the definitions of 
new coupling constants 

.9o = 9 , 9i = g' + 2cr , 

.92 = *■(/' -3cg' -3c 2 r)/ 5 , (2.30) 

and new control parameters 



f =T, fx = (r - r){g' - 2cr)/g, 

r 2 )/2g. (2.31) 



h = z+(r 2 



Note that the three control parameters fo, fx, and h go 
to zero in the mean-field theory of the collapse of large 
branched polymers, i.e., near z c , and that, however, the 
coupling constant g\ stays finite. 

The elimination of all terms that are at least irrelevant 
(in the sense of a naive scaling consideration) in compari- 
son to the retained ones reduces the stochastic functional 
to 



J' z 



d d xdt\{ h 



{nfA-^ + tfo-V 2 ) 

50 ~2 , 91 ~ 2 , 52 3" 

— n n H nn H n 

2 2 6 - 



r i 2 



hn 



}• (2-32) 



At this place we remark that for g2 — the stochastic 
functional J' z has the same form as the stochastic func- 
tional for ordinary directed percolation with a quadratic 
rapidity-reversal symmetry breaking ~ f\. Thus, in this 
special case the collapse transition belongs to the DP- 
universality class. Hence, we want to show wether gi is a 
relevant perturbation or it goes to zero under renormal- 
ization. 

The modifications of the dynamical functional includ- 
ing the elimination of the irrelevant higher order cou- 
plings make it worthwhile to discuss the support of the 
functional integration over the fields h and n anew. Like 
before, the integration of n, the density of monomers, 
goes along the positive real axes. We write the third 
order terms of the integrand of J z as 



fo - 2 , 5i ~ 2 

— n n H nn 

2 2 



^n 3 
6 



5o [n - 



5i 
— 1 

25o 



52 

3 



91 w 



4ffo 



(2.33) 



This form shows that the support of integration of h must 
be so that h — g^njlg^ is imaginary if this variable be- 
comes large. The convergence of the functional integral 
then requires 



5o > , 



45o52 



3(5o5i f > 



(2.34) 



to make sense beyond perturbation theory. If these con- 
ditions arc violated, higher order contributions become 
relevant for stability reasons, and a first-order transition 



takes place. Next, we consider the quadratic terms of J' z . 
Writing these terms as 



\\% - V 2 ]n + y n 2 = (h - *Ln) [f - V 2 ], 



51 









(fx + — f ) 




n 


IA g J 


50 - 





n (2.35) 



1 

— n 
2 



we see that we have to require 
fo > , gin + (<7 <7i)fo > , (g gi) > (2.36) 

for stability of the saddle-point at h = n = (for 
h = h c = 0). If these conditions are violated, a phase 
transition takes place. 

Through the 5-fimction condition in Eq. (|2.7p which 
selects from all the generated clusters only the clusters 
with given mass N, the stochastic functional J' z has lost 
its typical causal structure, whereas the absorbing state 
condition still holds, J^[n,n = 0] = 0. To restore the 
causal structure, we perform the duality transformation 

n(t) = itp(-t) , h(t) = -iip(-t) . (2.37) 

These steps lead us from Eq. (|2.32[) to 



J = f d d xdt\{(p X~ 1 d t + (f - V 2 ) 



50 - 2 , 5l-2 52 -3 



T l ~2 

-ih(p\. (2.38) 



Evidently, this new functional obeys causality, J\(p = 
0, (p] = 0, but the absorbing state property is satisfied 
only if h = T\ = g^ = 0. If only g^ = holds, it defines 
RFT with added linear and quadratic symmetry breaking 
terms. In the following, two facts will play an important 
role: (i) RFT is equivalent to ordinary DP, and (ii) the 
added symmetry breaking terms do not alter the fixed 
point structure of RFT. 

As it stands, the J is a truly proper field theoretic 
stochastic functional in the sense of being a minimal 
model. Nevertheless, it still contains that much gener- 
ality that allows to study the statistics of swollen DBP 
and the CDBP transition. These will be the subjects of 
the following sections. 



III. A BRIEF VIEW ON SWOLLEN DIRECTED 
BRANCHED POLYMERS 

In the previous section we have learned inter alia that 
the shape of V(N) for large N is controlled by the im- 
mediate vicinity of the branch point z c and that corre- 
lation lengths are large when the parameter r defined in 
Eq. (|2.18[) is small. For studying large clusters based on 
our functional J , Eq. (|2.38[) . this means that we should 
focus on the limit of vanishing control parameters h and 
fo- The collapse transition of these large polymers oc- 
curs for vanishing fx. This transition will be treated in 
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the later sections. Here we are interested in the swollen 
phase of these polymers, and hence we consider a fx that 
is positive and finite, say of the order of 1. Then, fi can 
be eliminated through a simple redefinition of the fields 



9o, 



<p t (fi/2)-VV -> (f 1 /2) 1 /% 
(fi/2) J / 2 gi -> 51, (fi/2) 3 / 2 g 2 -» .92, or formally by 
setting fx = 2. 

Now, we count engineering dimensions. As usual, we 
employ the some inverse length scale fi which is conve- 
nient for this task. Recalling that d denotes the number 
of transversal dimensions only, we have 



ip - 




2)/2 j 


<P ~ 


,Xd+2)/2 
H 1 ) 


(3.1a) 








9o r 


„ ,,(6-d)/2 


(3.1b) 


9i ' 


V 2 


-d)/2 

1 


92 


_ -(d+2)/2 _ 


(3.1c) 



We see that the upper critical dimension d c is 6 and 
that g '■= go is the only relevant coupling constant be- 
low d c - Neglecting the irrelevant couplings, we obtain 
the response functional 



Jyl = / d 



d xdt ^if^dfip + X^p-^ -A<^ 2 j, (3.2) 



Sip 

that is known to describe the relaxational dynamics at 
the Yang-Lee singularity edge [3l|. Here, 



i d x{\{Vipf 



, -9 3 



i/i^l . (3.3) 



is the Hamiltonian of an Isin g-or der parameter in an 
imaginary field near criticality [321 ] with r = To- It was 
shown [33l [3^ that the entropic exponent of the proba- 
bility V(N) ~ N~® exp(-Y/Y ) and the exponent of the 
transversal gyration radius R± ~ N"-*- 1 - are determined 
(as conjectured by Day and Lubensky [35[) by 



9=1 + (JYL ■ 

va.x = 0/ d > 



(3.4) 
(3.5) 



where (Jyl is the Yang-Lee exponent which relates the 
order parameter and the imaginary magnetic field via 



($-*c 



h — h. 



The (e = 6 — d)-expansion of 



the longitudinal Flory exponent is [3ll . [3 
1 

V A,\\ = r 



1 



95 
108 



9 4 

2 ln 3ie 



24+^ 



(3.6) 



Table U compiles values for 9 and v^y based on exact 
results for d = 0, 1 and d = 2 [3(| and the third order e- 
expansion results of de Alcantara Bonfim et al. [53] ■ For 
corrections to scaling and the derivation of the correction 
exponent in e-expansion see [HI, [39j . 



IV. FIELD THEORY OF COLLAPSING 
DIRECTED BRANCHED POLYMERS 

Now, we return to the CDBP transition as the main 
topic of our paper. 



d : 





1 


2 


3 


4 


5 


6 


: 





1/2 


5/6 


1.08 


1.26 


1.40 


3/2 


i>A,\\ ■ 


1 


0.79 


0.69 


0.63 


0.58 


0.54 


1/2 



TABLE I: Values for the critical exponents 8 and Vaa\ f° r 
various transversal dimensions d. 



Renormalization 



Our response functional J , Eq. (|2.38|) , describes specif- 
ically the CDBP transition when both control parameters 
fo and f\ become critical. In this case, they scale in terms 
of the inverse length scale [i as fo ~ fx ~ A* 2 - Then the 
engineering dimensions of the fields and the remaining 
parameters follow as 



V - / /2 , 
fo ~ f 1 ~ 
9o ~ 9i ~ .92 



h ~ /.t 

M (4-d)/2 _ 



(4+d)/2 



(4.1) 
(4.2) 
(4.3) 



Now, the upper critical dimension is d c = 4. All three 
coupling constants are relevant below d c . There are no 
further contributions that are relevant under the condi- 
tion of causality (originating from the absorbing condi- 
tion of the process) which we have restored through the 
duality transformation (|2.37p . 

It is well known that a perturbation expansion of cor- 
relation and response functions (generally called Greens 
functions) based on a field theoretic model like J pro- 
duces UV-singularities which have to be regularized and 
renormalized by singular counter terms. Here, we use 
minimal renormalization, i.e., dimensional regularization 
followed by minimal subtraction of e-poles, to cure the 
theory from UV-divergencies. For the general princi- 
ples and methods of renormalization theory see, e.g., 
Rcf. ESJ Ell- 
in the following, we use a ring ° to mark bare (un- 
renormalized) quantities, i.e., we let ip — > ip, and so on. 
From here on, quantities without a ring are understood 
as renormalized quantities. To get rid of £-poles, we use 
the renormalization scheme 



ip = Z 1/2 (ip + K<p) 
A = Z~ l Z x \, 



r = Z^g-f 
1 



h = Z 1 '^- 1 (h + ±G*/ V £/2 Z • A ■ t) 



9 a 



(4.4a) 
(4.4b) 

C-f, 
(4.4c) 
(4.4d) 



where e = 4 - d and G e = T(l + e/2)/(47r) d / 2 . Here, we 
have introduced the two-dimensional vector f = (fo,fi). 
Note that the mixing-term proportional to K reintro- 
duces the aforementioned gradient term, see the discus- 
sion below Eq. (|2.29|) . In minimal renormalization with 

dimensional expansion the additive contributions r , h c , 
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C_ become formally zero, and the e-content of the other 
covmtertcrms arc defined by pure Laurent-series: 



y(fc) 



(4.5) 



k=l 



l=k 



where the counterterms Y^' k ^ , when determined at the 
loop-order I, arc homogeneous polynomials of the renor- 
malized dimensionless coupling constants u a of order 21: 
Y^ l ' k \su a ) = s 2l Y^ l,k \u a ). Corresponding expansions 
are valid for the various other counterterms K, A, B a 
which themselves are homogeneous polynomials of order 
21, 21 — 1, 21 + 1, respectively, at each loop-order I. We 
calculate the various renormalization factors Z and ad- 
ditional counter-terms using a two-tiered approach, viz. 
we carry out an explicit 1-loop calculation and we employ 
Ward identities. The latter will be discussed further be- 
low. Our 1-loop calculation, see the Appendix for details, 
produces 



Z = 



Z = 



U Ul 



K 



Bn = 




Z\ = 1 



8e 



2ui —2uo 
3« 2 4ui 



u 2 , A 



uo I 1 
2e\l 0, 



u ui (2\u\ 



4u 2 — 3u 2 Uq 

4e 



u . 



, (4.6a) 
(4.6b) 

(4.6c) 

u u 2 , 
(4.6d) 

(4.6c) 



where ellipsis symbolize terms of higher loop-order which 
are regular expansions in the coupling constants u a but, 
of course singular in e. Note that in the case u 2 = 0, 
where rapidity inversion symmetry holds, K and B 2 are 
zero to all orders. In this case we have also Bq/uo = 
Bi/u\. In the general case all renormalizations are trivial 
if uo = because then only diagrams without loops are 
generated by the perturbation expansion. 



RG differential operator 

£>m = ^ + CA<9 A + £ ■ k ■ d± + E a & 9 «« 
( 7 /2 - Qh + Gl /2 n~ £/2 {t ■ a • r)/2 



d. 



(4.7) 



which defines the infinitesimal generator of the RG, the 
(half-) group of scale changes of the length The 
application of T>^ to renormalized fields in averages pro- 
duces 



7 - 



7 7 - 



(4.8) 



Recalling the renormalization scheme (|4.4|l and the form 
of the renormalization factors Z. and the additional 
counter-terms discussed below it, we find 



1 00 
M^lnZ..| = --(u- d u )YP = -£/y('.D , 

(4.9a) 

OO 
1=1 

l=l 

h( U .5 u )-i)sw = f;^ i ) ! 



= 2 



(4.9b) 
(4.9c) 
(4.9d) 



1=1 



where (u ■ d u ) = ^2 a u a du a , as well as 

C= (Ud^lnAlo =7~7a, 
t-k = M^llo = (7a -7) 



(4.10a) 
(4.10b) 



f3 a = fid^Ua^ = + ^ + 7 A ^) u a + b a , (4.10c) 



= ( 7 a - 1 /2)h + G 1 J 2 ^' 2 {t ■ a ■ £)/2 . 



(4.10d) 



for the RG functions. Explicit 1-loop results for these 
RG functions follow readily from Eqs. (|4.6p . 



C. Ward identities, redundancy, and the invariant 
renormalization group equations 



B. The renormalization group 

Next, we derive a renormalization group (RG) in the 
usual way by utilizing that the bare theory is independent 
of the external, arbitrary inverse length scale fi which 
enters through the renormalization scheme. Because of 
this independence, the derivative ^d/d^,Q,\o (|o denotes 
derivatives holding bare parameters fixed) vanishes for 
each unrenormalizcd quantity Q. Switching from bare to 
renormalized quantities Q depending only on renormal- 
ized parameters, the derivative /id/d[i\o changes to the 



It has been noted since the early days of field theory 
that each symmetry transformation (form-invariancc) 
which transforms not only fields but also parameters of a 
field theoretic functional implies the redundancy of one 
of its parameters and vice versa [42| . Redundancy means 
that a respective parameter is related to other parameters 
of the functional through simple linear transformations 
of the fields. These transformations, however, must not 
change the physical contents of the functional. Hence, in 
a proper field theory, redundant parameters can be and 
must be eliminated through linear transformations of the 
fields. 
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As we have already noted, the response functional J , 
Eq. (|2.38p . is form-invariant under the three continuous 
transformations (|2.28|) . The mixing and shift transfor- 
mations, Eq. (|2.29p , were used to eliminate two of the 
corresponding redundant parameters, r and c, and they 
introduced r = f o as a small free parameter. Now, we 
use this freedom to derive Ward identities from the shift 
transformation, Eq. (|2.28c|) . which reads 



(p — ► cp = cp + i<y , 



If -> Lp 



(4.11) 



when expressed in terms of the fields tp and (p. The re- 
sponse functional J, Eq. (|2.38| . is invariant under this 
shift, 

J^i^tiQaM = J{<p',tp';f,g a ,h'] , (4.12) 
if we also change the control parameters as follows: 

£^£'=£ + 7<7, (4.13a) 
h' = h- 1 r Q - ^-g . (4.13b) 



Here, we have introduced the two-dimensional vector 
9 = (<?0) — .9i )• Primarily, the shift is defined for un- 
renormalized quantities. However, because it does not 
involve a transformation of the coupling constants, and 
the renormalization constants are only functions of the 
couplings, the shift-invariance, Eq. (|4.12|) . also holds for 



-1/2 



U a fl 



e/2 



the renormalized quantities with g a := G t 

Now, wc compare the bare transformations with their 
renormalized counterpart. We renormalizc the bare form 

of Eq. (|4.13a[) . £ - £ = 75 with 7 = Z 1/2 j: 



i ■(£'-£) = 7G7 1 / V /2 (« + B) , 



(4.14) 



where u — (uq. — Ui) and B_ = (_Bq, — B\). Splitting this 
equation in singular and non-singular parts, we get 

£' =Z + jG- 1/2 n s/2 u, (4.15a) 
B=¥ = -u, (4.15b) 

with Y_ = Z_ — 1. The last equation, the first 
Ward-identity, is easily checked at 1-loop order by us- 
ing Eqs (|4.6[) . Wc perform the same procedure with 
Eq. (|4.13bj) and obtain 



Y 



-u- A. 



(4.16a) 
(4.16b) 



with Y_ = (lo,o,^o,i)- Also the second Ward-identity, 
Eq. (|4.16bj) . is easily checked at 1-loop order by using 
Eqs gU). 

The Ward identities (|4.15bj) and (|4.16bj) imply that 
there are further relations between the RG functions, viz. 







-(— e + j)u + u ■ k , 7 = w-ffl, 



(4.17) 



where 7 = (7001 7oi)- These relations can be used in- 
ter alia to check the consistency of the 1-loop results 
for the various RG functions and the results (|4.6|) lead- 
ing to them. Our results fulfill this consistency check. 
More generally we derive relations between different ver- 
tex functions. Note that the vertex function — Tj, , is 
diagrammatically defined as the sum of all one-line ir- 
reducible diagrams with k amputated external ip- and k 
amputated external (p-lines [4l| . Let T[ip, ip;r, g a , h] be 
the generating functional of the vertex functions, which 
we call the dynamic free-energy functional. It has the 
same invariance property, Eq. (|4.12[) . as the stochastic 
functional J ', which is identical to the mf-approximation 
of r. Hence, we have 

r[<P,<P;i>9u,h>] = r[<£, ¥?;£,#«, 0] + \h J d d xdt(p{-K,t) 
= ^[<P, ( P + iT,t + li>9a,h-jf -7 2 ,9o/2] , (4.18) 
Setting 7 = after differentiation leads to 



6(p{-x,t) V dr an dh> 



which is sometimes called equation of motion. This gen- 
eral property of the dynamic free-energy energy func- 
tional holds equally well in renormalized and unrenor- 
malizcd form, and is very helpful for higher-loop calcula- 
tions. Especially, we find by taking repeated functional 
derivatives 

r M (0) = fo + i(go J- - <7i J-)r li0 (0) , (4.20a) 



r 1 , a (o) = i( flb ^- fll A)r lll( o) ! 
r 2>1 (o) = i( 5 o^- 9l ^)r 2 , (o) ; 



(4.20b) 
(4.20c) 



where the argument means that the wavevectors and 
frequencies of the vertex functions are set to zero. 

Up to this point we have not used the scaling trans- 
formation, Eq. (|2.28ap . for eliminating a further redun- 
dancy. Applying this rescaling of the fields, J remains 
invariant if we change its parameters to 



TO — * T , 

.90 ->■ a~ 1 g , 



n 



an 



.91 ->■ "3i 



-> ah , 
.92 -> a 3 g 2 



(4.21) 



This invariance holds likewise for the renormalized and 
unrenormalized theory. Hence, it would be absolutely 
wrong to search at this place for fixed points, i.e., zeros 
of of all three Gell-Mann-Low functions, [5 a = 0, because 
the free scale transformation poisons the renormalization 
flow. Before we can search for fixed points, we first must 
eliminate the scaling redundancy through the definition 
of scale invariant quantities. We do that by defining the 
invariant dimcnsionless coupling constants 



U = UqUi , W — UqU2 , 



(4.22a) 
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and changing the fields and control parameters to their 
scaling invariant counterparts 



9a f , <P = 9of , 

,2, 



To = to , n = u h , h = g h. 



(4.22b) 
(4.22c) 



The RG differential operator X> M applied only to scaling 
invariant quantities reduces to 



fidfj, + CA9 A + t • k ■ dr + (3 u d u + Pwd n 



+ [hi/2-C)h + ( 2 



(4.23) 



and its application to the invariant fields in averages pro- 
duces 



T, , 70 , 72 7 

= -~2^~ ~2 < f > - 



2 V 



(4.24) 



Here, we have defined the invariant RG functions 

7o = 7 - e - Co , 7i = 7 + e + Co , 72 = "oV , ( 4 -25) 
where 



, 2/3 

Co — 

u 



(4.26) 



The invariant matrices k and a follow from their respec- 
tive counterparts k and a as 

K oo = koo , K io = u q 2 kw , 

Koi = ulk w , Ku = fcn + Co i (4.27a) 

"oo = u a 00 , a-io = «oi = lfl io j 

a n = Uo ~ 3 ai . (4.27b) 

From the above, we can straightforwardly collect the 
following explicit 1-loop for the invariant RG functions: 



7u 



7o 



7i 



72 = (j + ...)w, C 



5u 

T 

it 

~8 " 



Co 



and 



3m 

=" H h 

f 3u/8 3w/4 N 
^-1/2 19u/8-e y 

1/2 J 



(4.28a) 
(4.28b) 
(4.28c) 

(4.29a) 
(4.29b) 



In the DP case, iu = t\ = 0, the RGE has to have the 
usual form known from DP [18|, [2l[ . It follows that kqi 
and aoo have to be zero to all loop orders for w = 0. For 
the invariant Gcll-Mann-Low functions we find 



3 \ 3 

-e H — it u u; 

2 4 



37 

-2s + -, 



(4.30a) 
(4.30b) 



Once again, the Ward identities lead to relations be- 
tween the RG functions, namely 



too = (7- C) - ■ a) , «io = -(v-a)i, (4.31) 



and 



Pu/u 



(4.32) 



where we have introduced the two orthogonal vectors t; = 
(1, u _1 ) and u = (1, — u). We define 

Ko = (v ■ ^)o = «00 — UK10 = Kll - U~ Koi — flu/u, 

(4.33a) 

«1 = (^-«)o = «00 + K0lM _1 = Kll + K10U — [} U /U. 

(4.33b) 

Hence, if /3 U = 0, especially at a fixed point, it fol- 
lows from Eq. (|4. 32[) that v and v are a left and a right 
eigenvector of K, respectively, since the two vectors are 
orthogonal. In this case, we find the two eigenvalues 
Kq = —70/2 and K\. Of course, our 1-loop results for the 
invariant RG functions satisfy these relations holding to 
all loop-orders. 

To determine the critical behavior, we need to extract 
those quantities that are invariant under all the sym- 
metry transformations of the theory. Applying the shift 
transformation the dynamic free energy functional (|4.18D , 
that such a set of quantities is given by </> and the com- 
binations 



S = (f> + io , 

H = 2h - 2t(t - cr 2 , 

y = (r ■ v) = t + c , 



(4.34a) 
(4.34b) 
(4.34c) 



where a = T\ju = goTi/gi- Note that H is linearly 
related to the Laplace- variable z. 

In the following we are mainly interested in the average 
(S) = i(go^ + c) =: —iM since this quantity is linearly 
related to $ = (?i)dp- We recall that the critical part 
of $ yields the asymptotic behavior of the probability 
V(N) by an inverse Laplace transformation. Applying 
Eq. (|4.24[) together with (<j>) = 0, and using the relations 
(|4.3ip to (|4.33p between the RG-functions, we obtain the 
RGEs 



V^M = k q M 



_1 Koiy ; 



V„H = (n a + ~f X )H + (a 00 



2u 1 noi)y 2 ■ 



(4.35a) 
(4.35b) 
(4.35c) 



D. Fixed points and critical exponents 

Next, we search for the stable fixed points of the in- 
variant RGE P~2"4]) with the RG generator Eq. P~2"3"]) . 
The fixed point equations /3 U * = /3 W * = possess the fol- 
lowing solutions: (i) the trivial fixed point u„ = to* = 



12 



which is unstable below four transversal dimensions, (ii) 
a fully stable fixed point 



H~z-Zc 



0, 



(4.36) 



with stability (Wegner-)exponents X u = e + . . . and 
X w = 13e/12 + . . . (which correspond to the two leading 
correction-to-scaling exponents), and (iii) a fixed point 
with one stable and one unstable direction 



t/,„ 



16e 
~37~" 



10* 



13 

y 



13 

7 

12 



(4.37) 

which lies on a separatrix between the region of attrac- 
tion of the fully stable fixed point and a region where the 
renormalization flow tends to infinity, signalling possibly 
a discontinuous collapse transition. However, this fixed 
point lies in the unstable region, and is thus not accessible 
by our model. We conclude, therefore, that the renormal- 
ization group flow reaches asymptotically the fully stable 
fixed point. The line w = in the two-dimensional space 
of the invariant coupling constants has to be a fixed line 
of the renormalization flow because this line corresponds 
to the RFT-modcl, Eq. (|2~38|) with g 2 = 0, with its ra- 
pidity reversal invariance. Thus, W* = holds to all 
loop-orders. We conclude that as long as the collapse 
transition is continuous, it belongs to the universality 
class of directed percolation in all dimensions. 

We would like to underscore that this result, which is 
the main result of our paper, holds to arbitrary order in 
perturbation theory. The pivotal point is that for w — > 0, 
which is the case at the stable fixed point, J reduces to 
RFT and thus, the universal behavior governed by this 
fixed point is that of the DP universality class. The ex- 
plicit Tloop calculation was necessary to show the sta- 
bility of this fixed point. Higher orders in perturbation 
theory modify the two stability exponents quantitatively, 
and thus affect corrections to scaling. However, they do 
not change the stability of the fixed point qualitatively 
since in e-expansion, e is virtually an infinitesimal quan- 
tity. This does not answer, as always for the e-expansion, 
the question about a possible dimension d < d c at which 
the stability breaks down. 

We conclude this subsection by collecting the critical 
exponents of the CDBP transition. The anomalous scal- 
ing dimensions of the invariant fields are defined by 



f] = 70* + £ , V = 7i* - £ ■ 



(4.38) 



At the stable fixed point they coincide because (b = 
(2/3 /m )* = (Pu/u)* = 0. Also 72* = for w* = 0. 
Thus, the coinciding exponents are the same as for DP, 
i] = fj. The two eigenvalues kq* = — 7o*/2 = (e — rf)/2 
and >ti* define the order parameter and the correlation 
length exponents of DP, respectively, 



p/v = (2 -/so.) = (d + V )/2, 
1/v = (2 — Ku) . 



(4.39a) 
(4.39b) 




FIG. 3: Phase diagram with (y, H) as independent variables. 
The regions I and III are the physical parameter regions for 
a DP-process. The region III is excluded for the parame- 
ters of the statistics of directed branched polymers. The axes 
marked "red." correspond to some redundant variable such 
as gi/go, a, or c. The solid dot marks the DP-transition 
point corresponding to the collapse of large DBP. The line 
labeled H c (y) < with y > 0, where the order parameter 
$(z) shows a branching point singularity, correspond to large 
swollen DBP. The shading indicates region IV, where the 
order parameter 5>(z) becomes a complex function. At the 
straight line H c (y) = with y < between regions III and 
IV, the order-parameter shows essential singularities. The 
lines become surfaces if a redundant variable is admitted. 



7a* defines the dynamical exponent 

z = 2 + C* = 2 + 7} - 7a* ■ (4.39c) 
As is customary, we furthermore introduce the exponent 
A = p+ (z -rj)v, (4.39d) 
which is related to the scaling behavior of H . 



E. Scaling form of the equation of state and 
correlation lengths 

Next, we determine the equation of state. To this end, 
we introduce the dimensionlcss quantities M = il///i 2 , 
H = H/[i 4 , and y = y/n 2 , and seek M as a function of 
H and y, 



M = F(H,y) 



(4.40) 



At the stable DP fixed point, this relation between the 
three quantities must have scaling form. Using the 
critical DP exponents, Eq. (|4.39[> . and exploiting that 
Kqi = Q'oo = at the fixed line w = which contains 
the stable fixed point, we obtain from the flow equations 
(I4.35|) the fixed point flow 





V 


(4.41a) 


!*<*>- 




(4.41b) 


'■i 6 ® = 


V 


(4.41c) 
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with the RG flow-parameter £ defined by n(t) = 
The scaling form of the equation of state follows now 
from the solutions of the flow equations, and M(£) = 
F(H(£), $(£))■ After elimination of the flow parameter £, 
and the redefinitions M — > M, etc., we get the equivalent 
scaling forms 



where 



M 



Hf /A f±(y/\H\ 1/a ) = y"f(H/y*) , (4.42) 



,1/A 



where the subscript ± indicates the sign of H. Note 
that this field theoretic result is fully consistent with the 
mean- field equation of state as stated in Eq. (|2.18l) . 

Written in its scaling forms (|4.42[) . the equation of state 
leads to the following observations, see Fig. ([3]): the re- 
gion of the phase diagram corresponding to DP is given 
by M > and H > 0, where the primary control pa- 
rameter y can take both signs (43[. Note, however, that 
the stability condition, Eq. (|2.36|) restricts the variable y 
in our polymer problem to y > 0, and y = at the col- 
lapse transition, coinciding with the DP transition point. 
The primary control parameter here is the external field 
H ~ z — z c . The DP transition point becomes now a 
tricritical point given by H c = and y c = 0. The control 
parameter y constitutes here the crossover variable which 
drives the critical behavior from the DP-universality class 
to the critical point of the dynamical Yang-Lee model 
that we have discussed in Sec. IIII1 Hence, for y > we 
expect a critical line H = H c (y) < with M < in a 
region of the phase diagram not accessible for DP itself. 

To discuss the crossover in some more detail, we note 
that in mean- field theory (where (3 — 1 and A = 2), the 
scaling functions defined by Eq. (|4.42[) are simply given 
by 



f(x) = VxTT- l 



f± ( x ) = \/ X 2 ±\ — X . 



(4.43a) 
(4.43b) 



These equations show the tricritical point at x = and 
x = co. The critical line is described by the branchpoint 
singularity at x = — 1 or x = 1, respectively. However, 
mean-field theory for the crossover is applicable only for 
d > 6, where both types of polymers, collapsing and 
swollen, have mean-field behavior. 

To determine the general form of the crossover, we as- 
sume as usual that f{x) has in either case a singularity 
at same negative value x» of x (44j. To be specific, we 
assume that 



f(x) = f (x) + (x-x*) e ~ 1 f 1 (x) 



(4.44) 



fo(x) is the analytic part of f(x) at x — x ^ . m particular 
/o(a;*) = — A < 0. fi(x) may havefurthcr singularities 
at x = x* but it is assumed that /i(;c*) = B exists. 9 
denotes the entropic exponent of the swollen polymers, cf. 
Sec. IIIII With these assumptions, we find in the vicinity 
of the singular curve H = H c {y) = x*y A that 



M w -Af + By 



(H - x*y L 



(4.45) 



= l+/3/A 



(4.46) 



is the entropic exponent of the collapsing polymers. To 
conclude our discussion of the crossover, we emphasize 
that the problem involves two distinct kinds of nontrivial 
critical behavior with different upper critical dimensions. 
Thus, one cannot simply extract via analytical continua- 
tion e-expansion results for fix) from the known results 
for the DP equation of state [43||. A glance at the two- 
loop part of the latter shows, in deed, non-summablc 
singularities at x = — 1. 

Now, we turn to the scaling form of the correlation 
lengths, i.e., the diverging length-scales in transversal 
and longitudinal direction. To this end, we consider 
the scaling behavior of the correlation functions G k z of 

k fields <f> and k fields 0. The RGE for these correla- 
tion functions follows readily from the general RGE, see 
Eqs. (|4T2"3"]) and (g^l). At the stable DP fixed point with 
7o* = 7i* = 7* = f], 72* = 0, this RGE reads 



V^ + Mk + k) G kti ({x,t},y,H) = 0. (4.47) 



Via solving RGE (|4.47p . we obtain the scaling forms 
G.- k ({x,t},y,H) 



|ff| (fc+feWA / ±;M ({r/^,t/e||},2;/|i7r /A ) 

y {k+ ~ m h,%({*/UM\\},H/y A ) 



(4- 



of the correlation functions, where /, . fc k and f k k are 
scaling functions with the subscript ± denoting again the 
sign of H. Note that we need H < 0, for H < outside 
the physical region of DP. This scaling form implies the 
following scaling forms for the correlation lengths in the 
transversal and longitudinal directions, 



€x = \z-Zc\ v/A f±[y/\z- Zc\ 



,1/A 



WA /n (y/\ 



,1/A 



(4.49a) 
(4.49b) 



with scaling functions f± and f\\ with an analogous 
crossover behavior as discussed for the equation of state. 



V. ASYMPTOTIC PROPERTIES OF 
COLLAPSING DIRECTED BRANCHED 
POLYMERS 



Our field theory of Sec IIVI provided us with a number 
of results, in particular the scaling form of the equation 
of state and the correlation lengths. Though certainly 
interesting in their own right as they stand, they are not 
yet in a form or language that is customary in polymer 
theory. Here, we translate our results into results for the 
animal numbers and the radii of gyration. 
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As discussed in Sec. [TT1 the number of different poly- 
mer configurations A(N) for large N is up to an nonuni- 
versal exponential factor proportional to V(N), and 
the relation of the latter to $>(z) via inverse Laplace led 
asymptotically to Eq. (|2.19[) . Now we recall that &{z) is 
linearly related to M, and revisit our result (|4.42j) for the 
equation of state which we rewrite as 



$(z)-$(z c )~\z-z c f /A f(y/\z-z c \ 



,1/A 



(5.1) 



d : 


1 


2 


3 


4 


6 : 




1.10825 


1.27 


1.40 


1.5 


(e- 


Exp.) 


1.074 


1.255 


1.397 


3/2 


0: 




0.39151 


0.459 


0.49 


0.5 


(e- 


Exp.) 


0.375 


0.444 


0.486 


1/2 


"II ■ 




0.69007 


0.584 


0.54 


0.5 


(e- 


Exp.) 


0.601 


0.573 


0.539 


1/2 


v x ■ 




0.451494 


0.337 


0.285 


0.25 


(e- 


Exp.) 


0.376 


0.327 


0.285 


1/4 



Inserting this into the inverse Laplace transforma- 
tion (|2.19p . we obtain for N S> 1 the following scaling 
form for the animal numbers: 



TABLE II: Values of the CDBP exponents for various 
transversal dimensions d as produced by our field theory and 
numerical simulations 1271 14511 ■ 



A{N,y)~N- e f{yN+) 



(5.2) 



obtain 



with the animal exponent 6 which is the same as the 
entropic exponent defined in Eq. (|4.46[) . and the crossover 
exponent 



= 1/A. 



(5.3) 



We see here at the instance of the animal numbers, 
that the transcription of the field theoretic scaling re- 
sult to the corresponding scaling result in the polymer 
language asymptotically for N 3> 1 via inverse Laplace 
transformation amounts to simply replacing \z — z c \ by 
N -1 , see also the second part of Eq. (|2.19[) . In this way, 
after the inverse Laplace transformation of the correla- 
tion functions, space- and time-coordinates, which orig- 
inally scale with the correlation lengths, now scale with 
the appropriate powers of N . Consequently, we obtain 
for the longitudinal and the transversal radii of gyration 
as counterparts of the correlation lenghts 



V A,\\ 



e 

12 



f 37_ 



2 . 



£ 

18 
24 

e 

32 



53 

V288 ' 144 ^ 



ln(4/3) 
ln(4/3) 



13 


55 


288 


~ 144 


43 


17 


288 


~ 144 



(5.6a) 
(5.6b) 



(5.6c) 
(5.6d) 



In table HH we compare these e-expansions to results of 
numerical simulations as listed, e.g., in the review-article 
of Hinrichsen [45[ or the recent book of Henkel, Hinrich- 
sen, and Liibeck (27[, which, of course, are obtained by 
simulations of the DP-transition, and not for the polymer 
problem itself. Down to d = 2, the agreement between 
field theory and simulations is remarkably good. 



VI. CONCLUDING REMARKS 



R\\(N,y)~N v *.lf\\(yN+), (5.4a) 
R±(N,y)~N v ^f x (yN*), (5.4b) 

with critical exponents 

va,\\ = zv/ A , (5.5a) 
ua,± = v/A, (5.5b) 

where z is once again the dynamical exponent of DP 
which should not de confused with the variable z intro- 
duced through the Laplace transformation. 

It remains to state the e-expansions results for the 
polymer exponents encountered in this section. Using 
the relations of these exponents to the DP-exponents, 
which are known to second order in e [lU, [2^, [2l[, we 



In summary, we have developed a model for directed 
branched polymers in the framework of dynamical field 
field theory. As far as mean-field theory was concerned, 
our model reproduced a number of well known results 
and, therefore, passed important consistency checks. The 
same holds true for the field theory of swollen directed 
branched polymers which we extracted from our general 
model. Our main focus laid on the collapse transition 
of directed branched polymers for which, to our knowl- 
edge, no field theory existed hitherto. We showed to ar- 
bitrary order in perturbation theory that this transition 
belongs to the DP universality class. Because we used 
e-expansion, we can be sure about the stability of the 
DP fixed-point only for those dimensions for which the 
e-expansion is valid. Given that it is known from earlier 
work that the collapse transition of directed branched 
polymers belongs to the DP universality class in 1 + 1- 
dimensions and that our present work shows that this is 
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<0 =C(q,t) 



t) + — (0 = G(q,t) t) — 



FIG. 4: Diagrammatic elements. 







(a) 



(b) 



(c) 



(d) 



(e) 



(f) 



(g) 



(h) 




(i) 



FIG. 5: 1-loop contributions to the superficially divergent 
vertex functions. 



also the case in the vicinity of the upper critical dimen- 
sion 4+1, we are led to expect that this holds true in any 
dimension. We calculated the scaling behavior of several 
quantities of potential experimental relevance such as gy- 
ration radii and the probability distribution for finding a 
polymer consisting of N monomers. 
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APPENDIX A: 1-LOOP CALCULATION 

In this appendix we provide some details of our 1-loop 
calculation. The diagrammatic elements entering this 
calculation are the Gaussian propagator 

G(q, t) = Q(t) exp [-A (q 2 + f ) t] , (Al) 



the correlator 



C(q,i) 



n/2 

1 2 +t 



exp [-A (q 2 + f ) \t\] 



(A2) 



and the three-leg vertices —iXgo, —iXgi and iXg 2 , see 
Fig.UJ As usual, we concentrate on those vertex functions 
that are superficially divergent in the ultraviolet. These 
areT^i, r2,0j ^1,2, ^2,1, ^3,0 andr^o- The diagrammatic 
contributions to these vertex functions to 1-loop order are 
collected in Fig. [5] 

To give an example for the calculation of the various 
1-loop contributions, let us consider the diagram 3(c). 



The mathematical expression for this diagram is 

1 f°° 
3(c) = - (— iXgo)(iXg 2 ) / dt exp [- (iuj + 2Af ) t] 
Jo 

1^ J d d k exp [-A(k 2 + (q + k) 2 )t] 



(27T) 



(A3) 



where the last term indicates that this diagram can also 
be drawn with the arrows on the propagators pointing 
to the right which corresponds to the first term with ui 
replaced by — u>. The momentum and time integrations 
are straightforward. They lead to 



3(c) = - 



Xg g 2 r(e/2 - 1) 



4 (47r) d / 2 
+ (lu — > —u). 

e-expansion finally produces 



IU1 

2X 



7TT + T + -T- 



1 l-e/2 



3(c) = -Xg Q g 2 



Si 

4 



(A4) 



(A5) 



where we have omitted finite contributions that are in- 
consequential for our quest. Note that the contributions 
proportional to lu cancel as it should be the case. Note 
also that this diagram contributes a term proportional 
to q 2 to r2.o- A similar situation occurs for the collapse 
transition of isotropic randomly branched polymers. This 
and the renormalization mandated by it was overlooked 
in previous work [T2I [l3T | and hence led to incorrect re- 
sults. 

The remaining diagrams in Fig. [5] can be calculated 
by similar means. To save space, we merely list here the 
so-obtained results for the as yet unrenormalized vertex 
functions: 



f 1,1 = iuj + X (tq 
= iuj + X (fo 



3(a) - 3(6) 



, — e/2 G £ ( _ / iu) 

Ar ■^9o\9ot 1 - 9 i\to + - 



4 



(A6a) 



-Afi - 3(c) - 3(d) 



= — Afl + Afr 



-e/2 



G, 



90s gin +52 To + 



f 1,2 = iXg - 3(e) = iXg - iXf 6,2 — ^ g 2 gi 
f 2,1 = i\gi - 3(f) - 3(g) 

= i^9\ + iXf ~ e/2 5o {g g 2 - g\] , 
f 3,0 = -iXg 2 - 3(h) = -iXg 2 + iXr^ e/2 ^ 



(A6b) 
(A6c) 



ri, = 



(A6d) 
(A6e) 

iXh — 3(i) = iXh — iXf e ^ 2 go T o T i • (A6f) 
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It can easily be checked that these results satisfy the 
Ward identities discussed in Sec. IIVCI Application of 
the renormalization scheme (|4.4j) leads to renormalizcd 
vertex functions 

r li0 = Z 1 / 2 f lj0 , r 1)1 = Zf 1)l! (A7a) 

r 2 ,o = Z (f 2 , + 2KT 1,1) , r 1)2 = zV 2 t lt2 , (A7b) 

r 2il = Z^ 2 (f 2)1 + 2AT\ 2 ) , (A7c) 

r 3 ,o = Z 3 ' 2 (f 3 ,o + 3Kf 2A + 3K 2 f 1>2 ) , (A7d) 

and then finally to the 1-loop results for the renor- 
malization factors and additive counter-terms stated in 
Eq. gU). 

APPENDIX B: ESSENTIAL SINGULARITIES 





i 


r / v 












\ 

\\ 






















j \ 






















X 




r 



FIG. 6: (Color online) (d + 1) -dimensional cut of a droplet 
configuration (green inner part) of the true vacuum (s = 0, 
s = —2\r\/g) in the sea of the false vacuum (s = 2\r\/g, 
s = 0) and its approximation for ro > 1/ \/\t\ by a (d + 1)- 
dimensional cone (green inner and red outer part). 



In this appendix, we briefly discuss essential singular- 
ities appearing in the generating function $(2) for y < 
while crossing the separation line between regions III and 
IV in the phase diagram, Fig. [3J These singularities pro- 
duced by instantons are intimately related to the well 
known so-called false- vacuum problem [46[ . They lead to 
the behavior 



N' 6 ' cxp [-b(\y\V*N) 



CI 



(Bl) 



of the distribution function V with Q < 1 and a nonuni- 
versal amplitude b. Our discussion follows that of Har- 
ris and Lubensky [T3|, and Lubensky and McKane 
which are based on Langers classical paper [III on the 
droplet model. 

We have shown in the main text that the universal crit- 
ical behavior of collapsing directed branched polymers is 
equivalent to that of DP. Being interested in essential sin- 
gularities, we are hence led here to studying the behavior 
of the cluster distribution function of finite DP-clusters 
in the active phase above the transition point (where also 
an infinite cluster in time and space is generated with fi- 
nite probability). As the basis for this study, we use the 
stochastic functional (|2.38[) in its usual symmetrized real 
form with go = 9i — <?, 9i — 0, and with f\ eliminated. 
We set fo = t and work with 



J= I d d xdtX^s[\- 1 d t + T-W 



9 ( 
2^ 



- hsj . 

(B2) 

We need to understand the behavior of functional inte- 
grals with weight exp(— 3\s, s\) for r = — |t| near the first 
order line h w 0. To this end, we seek for saddle points. 
If h > near zero, the true vacuum, i.e., the saddle point 
producing the lowest value of J , is clearly the mean-field 
solution (s = 2|t|/<7, s = 0), up to corrections propor- 
tional to h. However, if we cross the border line at h = 
and go to small h < 0, the mean field solution becomes 
a false vacuum. At first sight, the saddle point solution 
(s = 0, s = —2\r\/g) appears to be the true vacuum. 



However, the response field has to satisfy the boundary 
condition s — * when the space coordinate or time tend 
to infinity. Thus, only non-stationary inhomogeneous so- 
lutions of the saddle-point equations, i.e., instantons or 
droplets, are allowed, see Fig. ©. 

To search for droplet solutions, it is useful to switch 
to dimensionless quantities. We introduce new fields s = 
2\r\<p/g, s = 2\r\(p/g, h = Ar 2 k/g 1 and scale time and 
space by \\r\t = t', It^^x = x' . In the following we will 
drop the primes for notational simplicity. The stochastic 
functional becomes 



9 

A= J d d xdt{(ftp- H(<p,<p)}, 
H(<p, (p) = -V<pV(p + (p(l - (p + (p)(p + kip . 



(B3a) 

(B3b) 
(B3c) 



The functional integral over exp(— J7[s, §]) is definitely 
equal to one. However, if we add a source term qs(x ,to) 
to the exponent, the saddle-point equations 



VV+ (1 - <p + 2<p)<p + k, 



-dp = V (p - 



'1 



WW: 



(B4a) 
(B4b) 



generate a droplet of the true vacuum that has its tip 
at (x ,to) and faces backward in time, see Fig. ([6]). We 
are interested in the critical droplet configuration whose 
lateral extension diverges as r — r c ~ with |fc| — > 0. 
We approximate this configuration by a critical spher- 
ical symmetric cone, see Fig. ([6]). Hence, we look for 
solitary waves tp(p) and <p(p) with p = r + vt — r c as 
kink-solutions of the saddle-point equations (|B4[) which 
have their variations near p = 0. By neglecting k as well 
as l/\x\ in comparison to v, we obtain from Eqs. (|B4p 
the differential equations 

V(p' = tp" + (1 - p + 2tp)tp, (B5a) 
-vip' = <p" + (1 + Cp - 2p)(p , (B5b) 
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FIG. 7: (Color online) The orbit of the droplet's instanton- 
configuration with |ar| = const, between the false and the true 
vacuum (green inner orbit), and its approximation by kink- 
solutions (red outer orbit). 



where the stroke means differentiation with respect to 
p. The properties of the kink-solutions, see Fig. ([7]) , and 
the fluctuations around them are extensively discussed 
in the work of Alessandrini et al. [4!| on classical kinks 
and their quantization in RFT. Also, in this work the 
appearance of the critical velocities v = — oo (the base of 
the cone) and v = 2 (the envelope of the cone) is shown. 



Taking all things together, the action A, Eq. (|B3b|) . 
of the large droplet consists of the usual terms arising 
from the volume ~ r d+1 of the droplet times k(p, where 
(p takes its value in the true vacuum, namely —1, and a 
surface part ~ r d arising from the variations of the kink- 
solutions (in the limit v — > 2 this surface part comes only 
from the base of the cone). We get, therefore, 



where a and b are positive constants. Minimizing A with 
respect to r, we obtain r — r c ~ l/|fc|, and finally 



Jc 



~ d/2 A 



(B7) 



A 



-a\k\r 



d+l 



br a 



(B6) 



Considering the eigenvalues of the Gaussian fluctua- 
tions about the critical instanton, we have (i) one nega- 
tive eigenvalue coming from the uniform variations of the 
extension r of the critical droplet leading to the imagi- 
nary part of $, (ii) [d+l) eigenvalues coming from 
the translations in space and time of the droplet and 
their compensation by the pinning of the droplet at its 
source, and (iii) a series of positive eigenvalues coming 
from the Goldstone modes of the broken translation sym- 
metries. For fully spherical symmetric droplets, the cal- 
culation of the latter eigenvalues is straightforward [5(| ■ 
In the present problem, however, the droplets are conical 
(spherical symmetric in space but directed in time). For 
this configuration, the Goldstone-mode eigenvalues have 
not yet been calculated, as far as we know, and their 
calculation represents an interesting and challenging op- 
portunity for future work. 

Finally, we return to the cluster distribution function 
V, Eq. (|Bip . which follows via inverse Laplace transfor- 
mation as discussed extensively in the main text. For the 
exponent (, this leads to £ = d/(d + 1). To determine 
9' , however, we first need to know the Goldstone-mode 
eigenvalues touched upon in the preceding paragraph. 
Thus, we here have to leave the interesting problem of 
calculating 9' for future work. 
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